EXACT SOLUTION OF THE SIX- VERTEX MODEL WITH DOMAIN 
WALL BOUNDARY CONDITIONS. FERROELECTRIC PHASE 



PAVEL BLEHER AND KARL LIECHTY 

Abstract. This is a continuation of the paper [4 of Blchcr and Fokin, in which the large 
n asymptotics is obtained for the partition function Z n of the six-vertex model with domain 
wall boundary conditions in the disordered phase. In the present paper we obtain the large 
n asymptotics of Z n in the ferroelectric phase. We prove that for any e > 0, as n — > oo, 
Z n = CG"F n [1 + 0(e~ n )], and we find the exact values of the constants C,G and 
F. The proof is based on the large n asymptotics for the underlying discrete orthogonal 
polynomials and on the Toda equation for the tau-function. 



1. Introduction and formulation of the main result 

1.1. Definition of the model. The six-vertex model, or the model of two-dimensional ice, 
is stated on a square n x n lattice with arrows on edges. The arrows obey the rule that 
at every vertex there are two arrows pointing in and two arrows pointing out. Such rule is 
sometimes called the ice-rule. There are only six possible configurations of arrows at each 
vertex, hence the name of the model, see Fig. 1. 

ii v v 



(1) (2) (3) 



(4) (5) (6) 

Figure 1. The six arrow configurations allowed at a vertex. 

We will consider the domain wall boundary conditions (DWBC), in which the arrows on 
the upper and lower boundaries point in the square, and the ones on the left and right 
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boundaries point out. One possible configuration with DWBC on the 4x4 lattice is shown 
on Fig. 2. 



Figure 2. An example of 4 x 4 configuration with DWBC. 

For each possible vertex state we assign a weight Wi, i = 1, . . . , 6, and define, as usual, 
the partition function as a sum over all possible arrow configurations of the product of the 
vertex weights, 

Zn= J2 W ^ w(*)=U W W=il V >? i{a) > 

arrow configurations a xdV n i=l 

where V n is the n x n set of vertices, t(x; a) e {1, . . . , 6} is the type of configuration a at 
vertex x according to Fig. 1, and Ni(a) is the number of vertices of type % in the configuration 
a. The sum is taken over all possible configurations obeying the given boundary condition. 
The Gibbs measure is defined then as 

M») = ^- (1-2) 

Our main goal is to obtain the large n asymptotics of the partition function Z n . 

The six-vertex model has six parameters: the weights By using some conservation 
laws it can be reduced to only two parameters. It is convenient to derive the conservation 
laws from the height function. 

1.2. Height function. Consider the dual lattice, 

V = {x = {i+ l -,j + \), 0<hj<n}. (1.3) 

Given a configuration cr on E, an integer-valued function h = h a on V is called a height 
function of a, if for any two neighboring points x, y G V, \x — y\ — 1, we have that 

%) - M^) = (1-4) 
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where s = if the arrow a e on the edge e G E, crossing the segment [x,y], is oriented in 
such a way that it points from left to right with respect to the vector xy , and s — 1 if a e is 
oriented from right to left with respect to xy. The ice-rule ensures that the height function 
h = h a exists for any configuration a. It is unique up to addition of a constant. Fig. 3 
shows a 5 x 5 configuration with a height function, and the corresponding alternating sign 
matrix, which is obtained from the configuration by replacing the vertex (5) of Fig. 1 by 1, 
the vertex (6) by (—1), and all the other vertices by 0. 
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Figure 3. A 5 x 5 configuration with a height function and the corresponding 
alternating sign matrix. 

Observe that if hi(x), h 2 (x), h 3 (x), h^{x) are the four values of the height function around 
a vertex x = (j,k), enumerated in the positive direction around x starting from the first 
quadrant, then the value of the element of the ASM is equal to 

hi(x) - h 2 {x) + h 3 (x) - h A (x) 
a jk = . (1.5) 



1.3. Conservation laws. Conservation laws are obtained in the paper [12] of Ferrari and 
Spohn, as a corollary of a path representation of the six-vertex model. We will derive them 
from the height function representation. Consider the height function h = h a on a diagonal 
sequence of points defined by the formula, 



Xj 

where both x and x^ lie on the boundary B' of the dual lattice V, 

S / = {x = (z + -,i), 0<i<n}U{x = {m+^,j + ~), < j < n} 

U{£ = (i + ~,n+~), 0<i<n}U{x=(-,j + -), 0<j<n}. 

Then it follows from the definition of the height function, that 

2, if t(x; a) = 3, 
h(xj) - h(xj-i) = I - 2, if t(x\ a) = 4 . 



where 



x 



0, if t(x; a) 



1,2,5,6, 



1.6) 



;i-7) 
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Hence 

= h(x k ) - h(x ) = 2N 3 (a; L) - 2N 4 (a; L), (1.10) 

where Ni(a; L) is the number of vertex states of type % in a on the line 

L = {x = x + (t,t), t e K}. (1.11) 

The line L is parallel to the diagonal y = x. By summing up over all possible lines L, we 
obtain that 

N 3 {a) - N 4 (a) = 0, (1.12) 

where Ni(a) is the total number of vertex states of the type % in the configuration a. 
Similarly by considering lines L parallel to the diagonal y = —x, we obtain that 

N x (a) - N 2 (a) = 0. (1.13) 

Also, 

N 5 (a)-N 6 (<r)=n, (1.14) 

which follows if we consider lines L parallel to the x-axis. 

The conservation laws allow to reduce the weights w±, . . . , wq to 3 parameters. Namely, 
we have that 

w^wpw^w^wpw? = C(n)a N ^a N *b N *b N *c N Sc N \ (1.15) 

where 

a = s/uhw^, b = y/w 3 w 4 , c = y/w 5 w 6 , (1-16) 

and the constant 

"M-sr <li7) 

This implies the relation between the partition functions, 

Z n (wi,w 2 , w 3 , w 4 , w 5 , w e ) = C(n)Z n (a, a, b, b, c, c), (1-18) 

and between the Gibbs measures, 

fi n (a; wi, w 2 , w 3 , w 4 , w 5 , w 6 ) = fj, n (a; a, a, b, b, c, c). (1-19) 

Therefore, for fixed boundary conditions, like DWBC, the general weights are reduced to the 
case when 

Wi = w 2 = a, w 3 = u> 4 = b, w 5 = w 6 = c. (1-20) 

Furthermore, 

Z n (a, a, b, b, c, c) = c n2 Z n (-, -, -, -, 1, l) (1.21) 

\c c c c J 

and 

f a a b b \ 

H n {a; a, a, 6, 6, c, c) = fi n I a; -, -, -, -, 1, 1 I , (1.22) 

V C C C C J 

so that a general weight reduces to the two parameters, -, - . 
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1.4. Exact solution of the six-vertex model for a finite n. Introduce the parameter 

a 2 + b 2 - c 2 

A = + - . (1.23) 

2ab v ; 

There are three physical phases in the six- vertex model: the ferroelectric phase, A > 1; the 
anti- ferroelectric phase, A < — 1; and, the disordered phase, — 1 < A < 1. In the three 
phases we parametrize the weights in the standard way: for the ferroelectric phase, 

a = sinh(t-7), 6 = sinh(t + 7), c = sinh(2|7|), < \>y\ < t, (1.24) 

for the anti-ferroelectric phase, 

a = sinh(7 — t), b = sinh(7 + t), c = sinh(27), \t\ < 7, (1-25) 

and for the disordered phase 

a = sin(7 — t), b = sin(7 + t), c = sin(27), \t\ < 7. (1-26) 

The phase diagram of the six-vertex model is shown on Fig. 4. 




1 a/c 



Figure 4. The phase diagram of the model, where F, AF and D mark ferro- 
electric, antiferroelectric, and disordered phases, respectively. The circular arc 
corresponds to the so-called "free fermion" line, when A = 0, and the three 
dots correspond to 1-, 2-, and 3-enumeration of alternating sign matrices. 



The phase diagram and the Bethe-Ansatz solution of the six-vertex model for periodic and 
anti-periodic boundary conditions are thoroughly discussed in the works of Lieb [21j-[24j, 
Lieb, Wu (25], Sutherland [30J, Baxter [2], Batchelor, Baxter, O'Rourke, Yung [3]. See also 
the work of Wu, Lin [31] , in which the Pfaffian solution for the six- vertex model with periodic 
boundary conditions is obtained on the free fermion line, A = 0. 

As concerns the six- vertex model with DWBC, it is noticed by Kuperberg [20J, that on 
the diagonal, 

a b , 

- = - = x, (1.27) 

c c 
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the six-vertex model with DWBC is equivalent to the s-enumeration of alternating sign 
matrices (ASM), in which the weight of each such matrix is equal to s N ~ , where AL is the 
number of (— l)'s in the matrix and s = . The exact solution for a finite n is known for 1-, 
2-, and 3-enumerations of ASMs, see the works by Kuperberg [20] and Colomo-Pronko [9] 
for a solution based on the Izergin-Korepin formula. A fascinating story of the discovery of 
the ASM formula is presented in the book [7J of Bressoud. On the free fermion line, 7 = f , 
the partition function of the six- vertex model with DWBC has a very simple form: Z n = 1. 
For a nice short proof of this formula see the work [9] of Colomo-Pronko. 

Here we will discuss the ferroelectric phase, and we will use parametrization fjl .24(1 . With- 
out loss of generality we may assume that 

7 > 0, (1.28) 

which corresponds to the region, 

b>a + c. (1.29) 

The parameter A in the ferroelectric phase reduces to 

A = cosh(2 7 ). (1.30) 

The six- vertex model with DWBC was introduced by Korepin in [16], who derived an 
important recursion relation for the partition function of the model. This lead to a beautiful 
determinantal formula of Izergin [13j for the partition function with DWBC. A detailed proof 
of this formula and its generalizations are given in the paper of Izergin, Coker, and Korepin 
[H]. When the weights are parameterized according to (11.241) . the formula of Izergin is 

_ [sinh(t - 7) sinh(t + 7 )]" 2 

(us*) 

where r n is the Hankel determinant, 
and 

m(A sinh(27) n 

smh(t + 7) smh(r — 7) 

An elegant derivation of the Isergin determinantal formula from the Yang-Baxter equation 
is given in the papers of Korepin, Zinn- Justin [19] and Kuperberg (20] (see also the book of 
Bressoud [7]). 

One of the applications of the determinantal formula is that it implies that the partition 
function r n solves the Toda equation 

2 d 
T n-r'n - r' n = r n+ iT n _!, n > 1, (') = — , (1.34) 



cf. the work of Sogo, [27]. The Toda equation was used by Korepin and Zinn- Justin [19] to 
derive the free energy of the six-vertex model with DWBC, assuming some Ansatz on the 
behavior of subdominant terms in the large N asymptotics of the free energy. 

Another application of the Izergin determinantal formula is that rjy can be expressed 
in terms of a partition function of a random matrix model and also in terms of related 
orthogonal polynomials, see the paper [32] of Zinn- Justin. In the ferroelectric phase the 
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expression in terms of orthogonal polynomials can be obtained as follows. For the evaluation 
of the Hankel determinant, let us write <f>(t) in the form of the Laplace transform of a discrete 
measure, 



m = . u , = 4Ve--sinh(2 7 /). 



sinh(t + 7) sinh(t — 7) ( 



1.35) 



Then 



where 



Tn = — J2 A(^) 2 n[ 2e ^ sinh (27^)], (1-36) 



h,...,l n =l »=i 



l<i<j<n 

is the Vandermonde determinant. 

Introduce now discrete monic polynomials Pj (x) = x J + . . . orthogonal on the set N = 
{/ = 1, 2, . . .} with respect to the weight, 

w(l) = 2e~ 2tl sinh(2 7 Z) = e ~ 2tw - e ~ m ~^\ (1.38) 

so that 

00 

J2 P j(l)Pk(l)™(l)=hkS jk . (1-39) 
1=1 

Then it follows from (11.361) that 

2 n2 l[h k , (1.40) 



n-l 



see Appendix in the end of the paper. We will prove the following asymptotics of hk 



Theorem 1.1. For any e > 0, as k — > oo ; 

^=(^Pr(i + o(^'"')). 

where 

q = e 2 ^ 2t . (1.42) 
The error term in {l.^l ) is uniform on any compact subset of the set 

{(t, 7 ): 0< 7 <t}. (1-43) 

1.5. Main result: Asymptotics of the partition function. This work is a continuation 
of the work [4] of the first author with Vladimir Fokin. In [4] the authors obtain the large n 
asymptotics of the partition function Z n in the disordered phase. They prove the conjecture 
of Paul Zinn- Justin [32] that the large n asymptotics of Z n in the disordered phase has the 
following form: for some e > 0, 

Z n = Cn K F n \\ + 0{n-% (1.44) 
and they find the exact value of the exponent k, 

1 2 7 2 



12 3n(n - 27) 



1.45) 
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The value of F in the disordered phase is given by the formula, 

_ 7r[sin(7 + t) sin(7 - t)] 
~ 2 7 cosg ' 

the exact value of constant C > is not yet known. 

Our main result in the present paper is the following theorem. 

Theorem 1.2. In the ferroelectric phase with t > 7 > 0, for any e > 0, as n — > 00, 



1.46) 



CG n F r 



1 + e 



1.47) 



where C = 1 — e 47 , G = e 7 *, and F = sinh(t + 7). The error term in (fi.^i| ) zs uniform 
on any compact subset of the set 

Up to a constant factor this result will follow from Theorem 11.11 To find the constant 
factor C we will use the Toda equation, combined with the asymptotics of C as t — > 00. The 
proof of Theorems 11.11 and 11.21 will be given below in Sections [2]l6j Here we would like to 
make some remarks concerning the phase transition between the ferroelectric and disordered 
phases. 

1.6. Order of the phase transition between the ferroelectric and disordered phases. 

We would like to compare the free energy in the disordered phase and in the ferroelectric 
phase, when we approach a point of phase transition. Consider first the ferroelectric phase. 
Observe that t, 7 — > as we approach the line of phase transition, 

b a 
- = - + 1, 



hence a, b, c — > in parametrization (11.241) . Consider the regime, 



In this regime, 

lim - = 

7^0 c 



t, 7-^+0, 
sinh(t + 7) a + 1 



t 

> a > 1. 

7 



lim 

7^0 sinh(27) 



lim — 

7^0 c 



sinh(t — 7) a — 1 



lim 

7^0 sinh(27) 



We have to rescale formula (11.471) according to (11.211) . 

'a a b b 



,1,1 

c c c c 



c~ n Z n (a, a, 6, 6, c, c) = OG n F n " 



1 + e- n 



where 



F sinh(t + 7) 

c 



Similarly, in the disordered phase, 



a a b b 
Z n ( ,1,1 

c c c c 



sinh(27) 
: Cn K Ff[l 



0(n~% 



where 



F Ti sin(7 — t) sin(7 + t) 

c 



27 sin(27) cos 



1.48) 
1.49) 

1.50) 

1.51) 
1.52) 

1.53) 
1.54) 



EXACT SOLUTION OF THE SIX- VERTEX MODEL 9 

Observe that parametrization (11.261) in the disordered phase is not convenient as we approach 
critical line (jl.48p . Namely, it corresponds to the limit when 

t, 7 ->--0, |— ->«>1. (1.55) 

1 2 7 

Therefore, we replace t for | — t and 7 for | — 7. This gives the parametrization, 

a = sin(t — 7), 6 = sin(t + 7 ), c = sin(2|7|), | 7 | < t. (1.56) 

The approach to critical line (11.481) is described by regime (11.491) . Formula (11.541) reads in 
the new t, 7 as 

7r sin(t — 7) sm(t + 7) 



We consider Fn on the line 



and we use the parameter 



on this line. In variables a, (3, 



(n — 27) sin(27) cos 
a + b 



2(f-7) 



(1.57) 



a, (1.58) 



b — a , 
(3 = (1.59) 



and 



F = — lYl ^ ne ferroelectric phase, (1.60) 
F = - — ^JL ' ^ in the disordered phase, (1-61) 



where 

Q(t,i) = - 

[71 — 27) sin 

forward calculation shows tha 

(3^1-0 



7r sin(t — 7) . « 

g(t, 7) = 1 A n (L62) 

_(tt-2 7 )_ 

A straightforward calculation shows that on the line = a in the disordered phase, as 



2(a - 1) 3 / 2 (1 - 8) 3 / 2 

git, 7) = 1 + 1 o J ^ W f + o((i - /?) 2 ). (1-63) 

37r(a; + l) 1 ^ 

By (11.601) . g(t, 7) = 1 in the ferroelectric phase. This implies that the free energy Fq exhibits 
a phase transition of the order | with respect to the parameter f3 at the point (3 — 1. Fig. 5 
depicts the graph of F = F Q (0) (the left graph) and its derivative, F^(f3) (the right graph), 
as a function of (3 = ^— on the line = 2. Observe the square root singularities of Fq at 
(3 = ±1, which correspond to the phase transition of order | . Since 

it is a phase transition of the order | as well, with respect to the parameter A at the point 
A = l. 
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The set-up for the remainder of the article is the following. In Section [2] we will discuss the 
Meixner polynomials, which will serve as a good approximation to the polynomials P n (z). In 
Section [3] we will discuss the Riemann-Hilbert approach to discrete orthogonal polynomials, 
and we will derive a basic identity, which will be used in the proof of Theorem ll.il In Section 
H] we will prove Theorem 11.11 Then, in Sections [5] and [6] we will obtain an explicit formula 
for the constant factor C, and we will finish the proof of Theorem 11.21 Finally, in Section [7] 
we will compare the asymptotics of the free energy in the ferroelectric phase with the energy 
of the ground state configuration. 



2. Meixner polynomials 

We will use the two weights: the weight w(l) defined in (11.381) and the exponential weight 
on N, 

W Q(Z) = q \ I g N; q = e^~ 2t < 1, (2.1) 

which can be viewed as an approximation to w(l) for large /. The orthogonal polynomials 
with the weight w®(l) are expressed in terms of the Meixner polynomials with (3=1, which 
are defined by the formula, 



M k ( Z;q )= 2 F l ^\- z ] i- q -^yj2 { -^T ] 



j=0 v ;j 

i -T j )'ncd(^-onS(*-<) 



(j'O 2 

j=o u ; 



Q 

i=o H 



see, e.g. [18]. For the corresponding monic polynomials, 



(2.2) 



They satisfy the orthogonality condition, 

oo —k£ 

J2M,(l;q)M k (l;q)q l = q -—^, (2.3) 



(2.4) 
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(M in Pjf stands for Meixner), the orthogonality condition reads 

oo / 7IS2 k 

£ Pf{l)P^W = hZ6 jk , hf = , ■ (2.5) 

1=0 ^ ^' 

They satisfy the three term recurrence relation, 

zPfiz) = P^(z) + **±h±l p?[ z ) +#±-p£ 1 {z\ (2.6) 



see [18]. According to (13. 8p . we take q = e 2j 21 . 

For our purposes it is convenient to introduce a shifted Meixner polynomial, 

(—1 ) k k^a k 

Q k (z) = P™(z - 1) = {1 _ q)k M k (z -l;q), (2.7) 

which is a monic polynomial as well. Equation (12. 5p implies the orthogonality condition, 

°° (k1) 2 n k+1 



By analogy with (ll.40p . define 

n— 1 

n ~r,2 

"■, 

k=0 

From KZSJl and (EH]) we obtain that 



2" 2 rK- ( 2 - 9 ) 



fLn2„fc+l on 2 „(n+l)n/2 

2" TT r } Lzi = / TT(^) 2 • (2.10) 



fc=0 v 11 1J k=0 

By analogy with (ll.3ip . define also 



z q _ [sinh( 7 + t) sinh( 7 - t)f 2 ^ Q 
ff(*!) a 

fc=0 

Then from (12.101) we obtain that 

Z% = F n2 G n , (2.12) 

where 

2 sinh(t - 7) sinh(t + ^)q 1/2 _ 2 sinh(t - 7) sinh(t + 7)e 7 ~* __ . . . , . 

^ ~~ fTT^ ~~ 1 _ e 2 7 ~2t ~~ smh{t + j), {Z.U) 

and 

G = q 1/2 = e^K (2.14) 

Our goal will be to compare the normalizing constants for orthogonal polynomials with the 
weights w and w®. To this end let us discuss the Riemann-Hilbert approach to discrete 
orthogonal polynomials. 
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3. RlEMANN HlLBERT APPROACH: INTERPOLATION PROBLEM 

The Riemann-Hilbert approach to discrete orthogonal polynomials is based on the follow- 
ing Interpolation Problem (IP), which was introduced in the paper [6] of Borodin and Bo- 
yarchenko under the name of the discrete Riemann-Hilbert problem. See also the monograph 
[5] of Baik, Kriecherbauer, McLaughlin, and Miller, in which it is called the Interpolation 
Problem. Let w(l) > be a weight function on N (it can be a more general discrete set, as 
discussed in [5] and [5], but we will need N in our problem). 

Interpolation Problem. For a given k — 0, 1, . . ., find a 2 x 2 matrix-valued function 
Y(z] k) = (Yij(z; fe))i<i i3 -<2 with the following properties: 

(1) Analyticity: Y(z; k) is an analytic function of z for z G C \ N. 

(2) Residues at poles: At each node I G N, the elements Yu(z; k) and Y2\{z\ k) of the 
matrix Y(z; k) are analytic functions of z, and the elements Y\2{z\ k) and Y^z; k) 
have a simple pole with the residues, 

Res Y j2 (z; k) = w^Y^l; k), j = 1, 2. (3.1) 

z=l 

(3) Asymptotics at infinity: As z —>■ oo, Y(z; k) admits the asymptotic expansion, 

(3.2) 



where D(z, r) is a disk of radius r > centered at z G C and 

lim n = 0. (3.3) 

It is not difficult to see (see [6] and [5]) that under some conditions on w(l), the IP has a 
unique solution, which is 

Y(z-k)-( Pkiz) C(wP k )(z) \ 

where the Cauchy transformation C is defined by the formula, 

W)=E4l. ( 3 - 5 ) 
i=i 

and Pk(z) = z k + . . . are monic polynomials orthogonal with the weight w(l), so that 

oo 

Y J P J {l)Pk{l)w{l) = h j 5 jk . (3.6) 
i=i 

It follows from fjiO]) . that 

h k = [Fi] 2 i, (3.7) 

where [Vi]2i is the (21)-element of the matrix Y\ on the right in (13.21) . In what follows we 
will consider the solution Y(z; k) for the weight w(l), introduced in fll.38p . 
Let Y® be a solution to the IP with the exponential weight w®, 

Y*(z-k\-( Qk(z) C(wQQ k )(z) \ 

{ ' } ~ {(hl^Q^z) {hti)- l C{w*Q k -i){z)) ■ ^ 
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Consider the quotient matrix, 

X{z;k) = Y{z;k)[Y Q {z;k)]- 1 . (3.9) 

Observe that detY^^z; k) has no poles and it approaches 1 as z — ► oo outside of the disks 
D(l, rf), I = 1, 2, . . ., hence 

det Y Q (z; k) = 1. 

Also, 

X(z; k) —> I as z — > oo outside of the disks n), / = 1, 2, . 
This implies that the matrix X can be written as 

X(z; k) = I + C[(w Q -w)R], 

where 

{h^^P^Qk^iz) -P k {z)Q k (z) 
Xhk-ih^-'Pk-iWQk^iz) -(hk^Pk^Qkiz) 
From formula (13.71) and (13.121) we obtain that 



R{z) 



(3.10) 
(3.11) 

(3.12) 
(3.13) 



hh — h 



Q 



— W( 



(3.14) 



i=i 



We will use this identity to estimate \hk — h^\. Observe that formula (13.121) can be further 
used to evaluate the large n asymptotics of the orthogonal polynomials P n (z), but we will 
not pursue it here. 

We would like to remark that identity (13.141) can be also derived as follows. Observe that 
since P^ and Qk are monic polynomials, the difference, — Qk, is a polynomial of degree 
less than k, hence 



J2 P k(l)[Qk(l)-Pk(l)]w(l) = 0. (3.15) 
i=i 

By adding this to equation (I3.6P with j = k, we obtain that 

oo 

h k = Y,Pk(l)Qk(l)w(l). (3.16) 
i=i 

Similarly, from (12.81) we obtain that 

oo 

^? = X)n(0Q*(0«' Q (0. (3-17) 
1=1 

By subtracting the last two equations, we obtain identity (13.141) . 

4. Evaluation of the ratio hk/h® 

In this section we will prove Theorem 11.11 By applying the Cauchy-Schwarz inequality to 
identity (13.141) . we obtain that 



J2Pk(l) 2 \w(l)-w^l)\ 



i=i 



1/2 



Y,Qk(l?\w(l)-w^l)\ 



1=1 



1/2 



(4.1) 
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— - 1 



From pg| . 



< 



' oo 1 V 2 r 1 oo 

tq E p *(0 2 HO - ™ Q (0 1 to E ^(o a HO - wQ ( 



1/2 



'fe Z=l 



|™(Z) - w Q (/)| = e~ (2 * +27)/ < GMO, Z > 1; Co 



e 4 ^ - 1 



hence 



' L k 1=1 

where 

Thus, by AO) , 
where 

By (USD, 



.. oo ^ oo 

■q E p *(0 2 MO - ™ Q (0 1 < ^o-Q E ^(0 VQ 



Ai fc 1=1 



C h k 



<C (l+e k ), 



£k 



h Q 



- 1 



e\ < C {l + £ k )S k , 

1 oo 

Sk = -s^2Qk(!)*\v>(!)-w Q ( 

h k 1=1 



1 X 

5k = TqE^WVo, 
^fc Z=l 



% = e 



-2(t+7) 



Let us evaluate <5fc. 

We partition the sum in (14. 8p into two parts: 



and 



where 



z=i 



^ = iq E g*(0Vo, 



L = [A; A ], < A < 1 
Let us estimate first We have from (12.71) . (12.81) that 

Q h (i) (-i) fe (i - q yi\^ 



(^) 1/2 



,1/2 



By O, 



M k (l - 1; g) = 1 + (1 - g-^(Z - 1) + (1 - g 



M fc (Z -l;g). 



_^ 2 A;(A;-l)(Z-l)(Z-2) 
(2!) 2 



_ fc(fc-l)(fc-2)(Z-l)(Z-2)(/-3) 
1 (3!) 2 
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If I < k, then the latter sum consists of I nonzero terms. For I < L it is estimated as 

M k (l - 1; q) = 0(k L L L+1 ) = o(e Llnk+{L+1)lnL ), (4.14) 

hence 

{h Qy /2 

Due to our choice of L in ( 14.111) . this implies the estimate, 

Q k (l) _ 0(e i^ +2fc , lnfc)> (416) 



Since < q < 1 and < A < 1, the expression on the right is exponentially small as k — > oo. 
From (14. 9p we obtain now that 

S' k = o( e feln(?+4fcAlnfc ). (4.17) 
Since A < 1 and q < 1, we obtain that 

= 0(e- c « fc ), co = > 0. (4.18) 

Let us estimate S' k \ 

By (EE}, 

1 oo 

r?yEO*(0V = i, (4-19) 



' l k 1=1 



hence 



-.CO s \ L oo 

™fc ;=L+1 \ ^ / Ai fc Z=L+1 



fc(0V<[*) =e~^ L . (4.20) 



=L+1 

Thus, 

& n < e -4 7 (^-l)_ ( 4 21 ) 

Since < A < 1 is an arbitrary number, we obtain from (14.181) and (I4.2ip that for any rj > 0, 

5 k = (e- kl - v ) . (4.22) 



Let us return back to inequality (14. 6p . Consider two cases: (1) e k > 1 and (2) e k < 1. In 
the first case (14.61) implies that 

e k < 2C 8 k , (4.23) 
which is impossible, because of ( I4.22p . Hence e k < 1, in which case (14. 6p gives that 

e\ < 2C 5 k . (4.24) 

Estimate (I4.22p implies now that for any i] > 0, 

e k = O (e- kl ~ v ) , (4.25) 



so that as k — ► oo, 

h k = h${l + e k ), |4| =e k = (e- fe "") . (4.26) 
This proves Theorem 11.11 
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From ( 14.26P we obtain that for any rj > 0, 



Z n = Z2l[(l + e h ) = CZ2 



k=0 



1 + e 



where 



C = JJ(1 + e k ) > 0. 



(4.27) 



(4.28) 



fc=0 



Thus, we have proved the following result. 
Proposition 4.1. For any e > 0, as n — > oo ; 



1 + e 



Z n = CF n G r 
where C > 0, F = sinh(t + 7), and G = e 7- '. 

To finish the proof of Theorem 11.21 it remains to find the constant C . 

5. Evaluation of the constant factor 



(4.29) 



In the next two sections we will find the exact value of the constant C in formula f)4.29p . 
This will be done in two steps: first, with the help of the Toda equation, we will find the 
form of the dependence of C on t, and second, we will find the large t asymptotics of C. By 
combining these two steps, we will obtain the exact value of C. In this section we will carry 
out the first step of our program. 

By dividing the Toda equation, fll.34p . by r^, we obtain that 

J'2 



T T" 
'n'n 



71' 



d_ 

dt 



The left hand side can be written as 



r t" 



,12 



^ =(lnr n )". 



From fll.40p we obtain that 



hence equation (15.11) implies that 



From (ll.4ip we obtain that 



We have that 



Ah r 



l>n-l 



T n +1 



(lnr n )" 



An 2 q 



i2n+l) 



4/i K 



tn-1 



+ O (e' n 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 
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hence from (15.41) . (15.51) we obtain that 

(lnr n )" = n 2 [- ln(l - e 2 ^ 2t )]" + O (e^) . (5.7) 
By (11.311) this implies that 

sinh(t — 7) sinh(t + 7) 



(lnZ„)" = n 2 
Since 



In 



I _ e 27-2t 



+ Ole- nl ~ e ). (5i 



. sinh(t — 7) sinh(t + 7) , r . , , , . , „ . . 

ln x _ e 2 7 -2i = ln[smh(t + 7)] + (t - 7) - In 2, (5.9) 

we can simplify (15.81) to 

{lnZ n )" = n 2 [lnsinh(t + 7)]" + O (V nl ~ £ ) . (5.10) 

Observe that the error term in the last formula is uniform when t belongs to a compact set 
on (7, 00), hence by integrating it we obtain that 

lnZ n = n 2 ln sinh(t + 7) + c x t + c + O (V" 1 "") , (5.11) 

where Co, c\ do not depend on t. In general, Co, c\ depend on 7 and n. By substituting 
formula (I4.29P into the preceding equation, we obtain that 

ln C + 71(7 - t) = c x t + c + O (V nl ~ £ ) . (5.12) 

Denote 

d = c — 717, di = c\ + n. (5.13) 

Then equation (15.121) reads 

lnC = dit + d + o(e- nl ~ E ). (5.14) 

Observe that C = C(7, t) does not depend on n, while dj = dj(j, n) does not depend on t, 
j — 1,2. Take any < 7 < t x < t 2 - Then 

ln C( 7 , t 2 ) - lnC( 7 , h) = d^h - h) + O (e" 1 ^) . (5.15) 
From this formula we obtain that the limit, 

lim di(7,n) = (^1(7), (5.16) 

n^oo 

exists. This in turn implies that the limit, 

lim d 2 (j,n) = 4(7), (5.17) 

n—>oo 

exists. By taking the limit n —>■ 00 in (I5.14p . we obtain that 

lnC r = di(7)t + d (T)- (5-18) 
Thus we have proved the following result. 

Proposition 5.1. The constant C in asymptotic formula \4-29$ has the form 

C = e dl(7)i+do(7) . (5.19) 
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6. Explicit formula for C 

In this section we will find the exact value of C, and by doing this we will finish the proof 
of Theorem 11.21 Let us consider the following regime: 

7 is fixed, t — > oo, (6.1) 

and let us evaluate the asymptotics of C in this regime. By (I3.6P and (11.381) we have that 

oo oo —24+27 — 2t— 27 

"0 = E »(') = E ^ tM ~" - ^ 2 "- 2 '") = ~ • < 6 - 2 > 

1=1 1=1 
Similarly, by f[2~B . 

e -2t+2 7 



^0 = 7—^7- > ( 6 - 3 ) 



hence 

Let us evaluate e k 
By gD, 



hh 

"k 



^ = + 0(e~ 2t ), t^oo. (6.4) 

for k > 1. 

4<G)(l + e fe )4, 0> = ^-j;- (6-5) 
In the partition of 5^ as + 5' k ' in (14.91) . (I4.10p . let us choose 

L=[0 3 +t 2 / 3 }. (6.6) 
From fl4TT2|) . fl4~13l) we obtain that for / < L, 

< { k-l),2 k L L L+l = e 2 7 -2 t ( g ?) 

hence 

n n n k ~ l h L T, L+1 n kh,LjL+l 

ri ^ yoy ft , y -27-2* \ 

5fc < : < — : , yo = e 7 . (6.8) 

1 - 9o 1 - yo 

In addition, by (jPDjl . 

<f£ < e" 47L . (6.9) 

Our choice of L in (16.61) ensures that there exists t > such that for any t > t and any 
k > 1, 

4 = ^ + ^<e- fel/2 -* 1/2 . (6.10) 
From (16. 5p we obtain now that for A; > 1 and large t, 

e k < de-^'^ , C 1 = (2C ) 1/2 . (6.11) 

By (gjZBD, 

00 

lnC = ^ln(l + £ fe ), |e*|=e fc . (6.12) 

k=0 

From equations (16 .4p and (16. lip we obtain now that 

lnC = ln(l-e~ 47 ) + 0(e-^), t 00. (6.13) 
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On the other hand, by (15.14p 



This implies that 
so that 



InC = di(7)f + do(7) 
di(7) = 0, rf (7) = ln(l-e- 4 ^) 
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(6.14) 
(6.15) 
(6.16) 



(7 = 1- e" 4 7 

By substituting expression (I6.16P into formula (I4.29p . we prove Theorem 11.21 

Let us compare now the asymptotics of the free energy in the ferroelectric phase with the 
energy of the ground state. 

7. Ground state configuration of the ferroelectric phase 



— e — ; — a— 



Figure 6. A ground state configuration. 

The ground state is the configuration 

(T5 if x is on the diagonal, 

a 3 if x is above the diagonal, 

cr 4 if x is below the diagonal, 

see Fig. 5. The weight of the ground state configuration is 



a gs (x) 



2 (C 



where 



w{a gs ) = b n 



F = sinh(t + 7 ), G 



F n2 G n 



o- 



sinh(27) 
sinh(t + 7) 



(7.1) 

(7.2) 
(7.3) 
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The ratio Z n /w(a gs ) is evaluated as 

Z 



n 



Gi, (7.4) 



where 



Observe that 



G _ e^sinh(t + 7 ) _ - e~ 2t 
Gl "C7 " sinh2 7 "e^-e- 2 ^ {f ' b) 

hm^= hm lnw( ; gS) =lnF, (7.6) 

so that the free energy is determined by the free energy of the ground state configuration. 
This can be explained by the fact that low energy excited states are local perturbations of the 
ground state around the diagonal. Namely, it is impossible to create a new configuration by 
perturbing the ground state locally away of the diagonal: the conservation law Ns(a) = N^a) 
forbids such a configuration. Therefore, typical configurations of the six- vertex model in the 
ferroelectric phase are frozen outside of a relatively small neighborhood of the diagonal. 

This behavior of typical configurations in the ferroelectric phase is in a big contrast with 
the situation in the disordered and anti-ferroelectric phases. Extensive rigorous, theoret- 
ical and numerical studies, see, e.g., the works of Cohn, Elkies, Propp [8], Eloranta [XT] . 
Syljuasen, Zvonarev [28], Allison, Reshetikhin pQ, Kenyon, Okounkov [15], Kenyon, Ok- 
ounkov, Sheffield [17], Sheffield [26], Ferrari, Spohn [12], Colomo, Pronko [10], Zinn- Justin 
[33] , and references therein, show that in the disordered and anti-ferroelectric phases the 
"arctic circle" phenomenon persists, so that there are macroscopically big frozen and ran- 
dom domains in typical configurations, separated in the limit n — > oo by an "arctic curve 



J) 



Appendix A. Derivation of formula ( 11.401) 

Multilinearity of the determinant function, combined with the form of the Vandermonde 
matrix, allows us to replace A(/j) with 



det 



/ 1 



1 

P 2 (/ 2 ) 



1 

P2(h) 



1 

Pl(lr 
P2(!r. 



\ 



(a.i; 



Pn-l(ln)/ 



where {Pj{x)}^L Q is the system of monic polynomials orthogonal with respect to the weight 
w{l). Then (11.361) becomes 



Tr, 



2 n 
n! 



h,...,l n =l fc=l / fe=l 



(A.2) 
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Note that the orthogonality condition ensures that, after summing, only diagonal terms are 
non-zero, so we get 

r)n 2 oo / n \ n n— 1 

r « = ^T E Ell^-it 1 *) mtKfc) =2» a n*»- (A.3) 

«i,...,/„=l \7r6S„ fc=l / fe=l fc=0 
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